Model of correlated sequential adsorption of colloidal particles 



R. Pastor-Satorras^ and J. M. Rubi^ 
^Dept. de Fisica i Enginyeria Nuclear, Universitat Politecnica de Catalunya, 
Campus Nord, Mddul B4, 08034 Barcelona, Spain 
^ Dept. de Fisica Fonamental, Facultat de Fisica — CER Fisica Sistemes Complexos, 
Universitat de Barcelona, Av. Diagonal 647, 08028 Barcelona, Spain 
(February 1, 2008) 

We present results of a new model of sequential adsorption in which the adsorbing particles are 
correlated with the particles attached to the substrate. The strength of the correlations is measured 
by a tunable parameter a. The model interpolates between free ballistic adsorption in the limit o 
oo and a strongly correlated phase, appearing for cr ^ and characterized by the emergence of highly 
ordered structures. The phenomenon is manifested through the analysis of several magnitudes, 
as the jamming limit and the particle-particle correlation function. The effect of correlations in 
one dimension manifests in the increased tendency to particle chaining in the substrate. In two 
dimensions the correlations induce a percolation transition, in which a spanning cluster of connected 
particles appears at a certain critical value CTc. Our study could be applicable to more general 
situations in which the coupling between correlations and disorder is relevant, as for example, in the 
presence of strong interparticle interactions. 

PACS numbers: 68.43.-h, 05.70.Ln 



I. INTRODUCTION 

The study of the irreversible adsorption of colloidal 
particles onto a surface has attracted a great deal of in- 
terest in the last years, due to its many practical appli- 
cations in physics, chemistry, biophysics, medicine, etc. 

The understanding of these processes has been 
deepened mainly through the formulation of different 
models, defined via a set of rules by which the particles 
accommodate upon arriving at the surface. The different 
rules are responsible for the different values of the quanti- 
ties describing the geometry of the absorbed phase, as the 
maximum fraction of the surface covered by particles — 
the jamming limit Ooo — or the particle-particle pair cor- 
relation function g(r). In the random sequential adsorp- 
tion model (RSA) particles are placed at randomly 
selected positions on the surface. When an incoming par- 
ticles overlaps with a previously adsorbed one, it is re- 
jected; otherwise, it becomes irreversibly adsorbed. The 
RSA model is thus a good approximation when particles 
arrive at the surface purely by diffusion, and excluded 
volume effects are predominant j^] . In the ballistic model 
(BM) [p|-pT[|, on the other hand, particles descend to the 
surface following straight vertical trajectories. An in- 
coming particle that does not reach the surface directly 
is allowed to roll over the previously adsorbed ones, fol- 
lowing the steepest descent path, until it reaches a stable 
position. Particles that eventually rest on the surface are 
irreversibly adsorbed; otherwise they are rejected. The 
BM is thus a valid approximation to describe adsorption 
in the presence of very strong interactions between par- 
ticles and substrate |I2[-fi4[|. 

In the models described above, and in many of the 
variations of them analyzed so far, particles are supposed 
to interact only through short-range interactions — hard- 



core repulsion. The main exceptions are the analysis 
of the role played by electrostatic and dipolar p^ ] 
interactions. In the absence of long-range interactions, 
the particles arriving at the surface are essentially un- 
correlated from the adsorbed phase, and interact with 
them only via excluded-volume effects. Therefore, little 
is known about the general effect of correlations among 
the adsorbed phase and the adsorbing particles ||l^ . 

Our purpose in this paper is to analyze the influence 
of correlations among particles in the structure of the ad- 
sorbed phase in a simple numerical model. In our model, 
particles adsorb sequentially onto the surface. The pres- 
ence of correlations affects the position on the surface in 
which the particles try to adsorb. In order to mimic the 
effect of long-range interactions, which tend to attract the 
incoming particles in the vicinity of the particles already 
attached, the trial position of the next particle is selected 
to be the position of a particle in the substrate plus an 
increment ^, selected at random from a given probability 
distribution p(^). Once the trial position has been se- 
lected, the adsorption process proceeds according to the 
rules of the BM model. The distribution p(^) depends 
on a parameter a that interpolates between a flat distri- 
bution and a delta function, which allows to explore the 
range between absence of correlations and a very strong 
correlation effect, respectively. Apart from the obvious 
relation with the adsorption in presence of long-range at- 
tractive interactions |14| , the inclusion of general correla- 
tions in adsorption phenomena is relevant in other kind 
of related problems, such as car parking, bird nesting, or 
adsorption with memory effects, where the state of the 
adsorbed phase at a given time exerts a strong influence 
on the position of the next incoming unit. We want to 
note the our model differs from the generalized ballistic- 
deposition models presented in Refs. due to the 
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FIG. 1. Adsorption rule for the correlated adsorption 
model. An adsorbed particle j is randomly selected. The 
newly adsorbed particle is deposited at a distance ^ from the 
center of the selected particle. If it does not overlap with any 
other adsorbed particle, it is irreversibly attached. Otherwise, 
it is adsorbed following the prescription of the BM deposition. 

implicit role as nucleation centers played in our model by 
the adsorbed particles. 

The paper is organized as follows. In Section II we in- 
troduce the model in a general context, whereas in Sec. 
Ill we present the results of numerical simulations in di- 
mensions d = 1 and 2. Finally, in the last Section we 
discuss the common features induced by correlations in 
the structure and geometry of the adsorbed phase. 

II. MODEL 

Our model considers the sequential adsorption of parti- 
cles of diameter a onto a d-dimensional surface; it is thus 
d + 1 dimensional, including the direction perpendicular 
to the plane in which the particles move before adsorp- 
tion. Particles arrive at the surface following straight 
vertical trajectories, and upon touching the surface or 
any adsorbed particle, they are attached following the 
rules of the BM. The key ingredient in the model, that 
differs from previous works, is the selection of the landing 
point of the arriving particles. In order to introduce cor- 
relations among the incoming particles and the adsorbed 
phase, the landing point is selected as follows: The first 
particle is locate at random on the substrate. Suppose 
that, at a given time step, the substrate is composed by 
set of N particles, located at positions Hi, i = 1, ... TV. 
The landing position of the next particle is chosen as 

R' = Rj+nC, (1) 

where Rj is the position on the surface of a particle j 
randomly selected among the N particles present on the 
adsorbed phase, n is a unity vector, parallel to the ad- 
sorption plane and oriented at random, and the distance 
^ is Gaussian random variable, distributed according to 
the density 

P(0 = ^^expf~-^) , (2) 
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FIG. 2. Chain-length density function P{£) for different 
values of the parameter ct in 1 -I- 1 dimensional correlated 
adsorption. 

see Fig. 0. The only parameter in the model is a, defined 
as the root-mean-square distance from an adsorbed par- 
ticle at which the new particle is probed. Once selected 
the position R', the particle is attached following the 
rules of the BM. If the particle is accepted, the updated 
position R" becomes the new value RAr+i. If the particle 
is rejected, a new position is chosen and the procedure 
iterated. 

The above algorithm has a very intuitive interpreta- 
tion. Each one of the particles in the substrate play the 
role of a nucleation center, attracting incoming particles 
towards them. We can envision this process as the re- 
sult of long-range attractive interactions between the ele- 
ments of the adsorbed phase and the incoming particles. 
The attraction has a random component given by the 
Gaussian term in Eq. (|]), which represents the effect of 
the forces exerted on the incoming particle by the regions 
of the substrate located far away from the landing point. 
The strength of the attraction is given by the parameter 
a in Eq. . With this parameter we can interpolate be- 
tween absolute lack of correlations, achieved in the limit 
(T — > cx), and the case of very strong correlations, in the 
limit (7^0. In the former case, the Gaussian (||) tends 
to a flat distribution, which corresponds to a free BM 
adsorption. On the other hand, in the limit cr — > the 
distribution (|^) becomes a sharply peaked delta function. 
In this case the correlations are maximal, and one would 
expect to obtain a close-packed substrate, at least in 1 -|- 1 
dimensions. 

In the following Section we report the result obtained 
from simulations performed with the random increment 
distribution given by Eq. (H). Other sets of simulations, 
done with a different peaked distribution [such as an 
exponential p{£,) = exp(— ^/o(T)/acr], yield qualitatively 
similar results, confirming the fact that the important 
point is the presence of correlations, and not its particu- 
lar form. 
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FIG. 3. Jamming limit O^o as a function of the param- 
eter a in 1 + 1 dimensional correlated adsorption. The 
dashed line signals the jamming limit for free BM adsorption, 
61™ ~ 0.80865. 



III. RESULTS 

In order to check the effects of correlations in the ad- 
sorption process, we have performed extensive numerical 
simulations of our model in 1 -I- 1 and 2-1-1 dimensions. 
Simulations were done in systems of different size L. Sta- 
tistical averages were performed over at least 100 differ- 
ent realizations. In order to avoid cross-effects between 
the characteristic length aa and the system size L, we 
always consider aa L. 



A. Adsorption in 1 + 1 dimensions 

We have performed simulations of the correlated ad- 
sorption model in 1 + 1 dimensions over a substrate con- 
sisting of a line of length L = 5000a, with periodic bound- 
ary conditions. 

In 1 -)- 1 dimensions, the most relevant feature of the 
model is the increased tendency of the adsorbed parti- 
cles to form connected structures, identified with chains, 
whose average length tends to increase with decreasing 
a. We can use as a measure of this kind of order the 
chain-length density function, P{£), defined as the av- 
erage number of chains of length £a per unit length of 
substrate. In Fig. ^ we represent the chain-length den- 
sity function computed at different values of a. The value 
a = oo corresponds to simulations of free BM adsorption. 
For large a we observe a very fast (super-exponential) 
decay of P{£), indicative of a lack of any characteristic 
length scale. For small values of a, on the other hand, 
P{£) shows a clear exponential tail, P{£) ^ exp(— 
with a characteristic length £c depending on a. In par- 
ticular, for <T — 1 we estimate a value £c = 2.48 ± 0.05. 

Another quantity of interest is the jamming limit Ooo, 
which in the present case can be related to the chain- 
length density function through the expression 
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FIG. 4. Particle-particle correlation function g{r) for dif- 
ferent values of a in 1 -|- 1 dimensional correlated adsorption. 
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In Fig. ^ we represent the jamming limit as function of 
the parameter a. For values of a larger than about 8 
we recover with good accuracy the limit of free BM ad- 
sorption, 0^ ~ 0.80865 [|. For smaller values of a, on 
the other hand, the presence of the correlations induces 
a higher compaction on the substrate, with a jamming 
limit that approaches the close-packing for a value of a 
equal to 0.5 {9oo = 0.9918 ± 0.0005). 

As a final probe of the substrate's structure, we have 
measured the particle-particle pair correlation function 
g{r). In Fig. ^ we plot g{r) as a function of the reduced 
distance r/a corresponding to six different values of a. 
We observe that the correlations enhance the maxima 
of g{r), which occur at distances pa for integer p (cor- 
responding to interparticle distances equal to a multiple 
of the diameter). The decay of the maxima is close to 
exponential again. For large values of a we recover the 
behavior of the free BM adsorption, with peaks decaying 
super-exponentially. 



B. Adsorption in 2 + 1 dimensions 

Simulations in 2 + 1 dimensions were performed on 
squares of size L = 125a, with periodic boundary condi- 
tions. In Fig. ^ we plot the jamming limit 9oo obtained 
from simulations for different values of a. The trend 
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FIG. 5. Jamming limit 6*00 as a function of the param- 
eter a in 2 + 1 dimensional correlated adsorption. The 
dashed line signals the jamming limit for free BM adsorption, 
61™ ~ 0.610. 



previously observed in 1 + 1 dimensions is here recov- 
ered. For large values of a the jamming limit tends to 
the free BM adsorption 9^ ~ 0.610 fill. For decreasing 
a, 9ac increases approaching the close-packing regime; 
for a = 0.5, for example, we estimate 6*00 = 0.715 ±0.001 
p9[ . It should be noted that the convergence of 6*00 to- 
wards its maximum value is slower in this case than in 
the adsorption in 1 -f 1 dimensions (compare Figs. ^ and 
1). 

In Fig. H we plot the particle-particle correlation func- 
tion g{r) corresponding to the 2-1-1 case. For a = 00, 
we recover the expected form in free BM adsorption (the 
divergent peak at r = a, corresponding to the close con- 
tact of particles, has been omitted for the sake of clar- 
ity). The secondary peak at r = 2a is clearly visible, as 
well an intermediate, much smaller peak, located approx- 
imately at r ~ 1.73a. For decreasing a the height of the 
secondary peak increases notably, and the intermediate 
peak becomes more and more noticeable. 

From the analysis of the correlation function we can 
conclude that the correlations in the adsorption mecha- 
nism induce the creation of connected structures on the 
substrate, which are responsible for the enhancement of 
the peaks in g{r). These connected structures cannot be 
characterized in terms of chains, since in 2 -I- 1 dimen- 
sions particles adsorb in multiply connected configura- 
tions, and therefore the notion of chain loses its meaning. 
In this case, however, the structure of the adsorbed phase 
can be better analyzed in terms of its percolation prop- 
erties ||2^,|l^,^. Since the particles can roll one above 
other before adsorption, the substrate becomes eventu- 
ally composed by clusters of connected particles, which 
can be easily identified. In free BM adsorption, the satu- 
rated phase has low connectivity and it is below the per- 
colation threshold: in simulations on any finite system of 
size L, the largest cluster of connected adsorbed particles 
is not a spanning cluster. In order to investigate the per- 
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FIG. 6. Particle-particle correlation function g{r) for dif- 
ferent values of a in 2 -f 1 dimensional correlated adsorption. 
The divergent peak at r = a has been suppressed for clarity. 

eolation properties of the correlated adsorption model, 
we have computed the susceptibility x of the distribu- 
tion of clusters of connected particles as a function of a. 
The susceptibility is defined by p3] 



X 



(4) 



where Us denotes the density of clusters of size s and 
the prime in the summation indicates that the biggest 
(spanning) cluster is excluded in the sum. At the perco- 
lation threshold CTc, the susceptibility is expected to ex- 
hibit a peak, with power-law decay both below and above 
threshold. In Fig. we have plotted the susceptibility x 
computed as a function of the parameter a. The presence 
of a maximum in located approximately at ac — 2.4, 
signals the presence of a percolation transition. Below 
the threshold, the clusters of connected particles are nu- 
merous and rather small. Above the threshold (cr < CTc), 
however, there exists an spanning cluster, that crosses the 
system from boundary to boundary. We show the pres- 
ence of such a spanning cluster in Fig. |[ For a > 2.4, 
the largest cluster in a typical simulation of our model is 
rather small. For a — 2.4 a largest cluster develops that 
first spans the whole system. For even smaller values of 
cr, the largest cluster covers most of the substrate. 
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FIG. 7. Susceptibility for the 2 + 1 dimensional correlated 
adsorption. The peak, located at cTc — 2.4, signals the pres- 
ence of a percolation transition. 



IV. CONCLUSIONS 

In this paper we have presented a new model of cor- 
related sequential adsorption of particles on a substrate 
which takes into account the existence of certain degree 
of correlation between the final positions of the adsorb- 
ing particles and the location of the preadsorbcd ones. 
Its implementation is carried out by modifying the rules 
of the BM, by imposing that the final position of an in- 
coming particle is randomly distributed, with a specified 
probability density, around the location of a particle in 
the substrate. In this way, the model can describe many 
different physical situation ranging from the absence of 
correlations, corresponding to BM model, to the pure 
deterministic case in which the distribution of particles 
practically reaches the close-packing configuration. 

The main consequence of the existence of correlations 
is the appearance of long-range order in the adsorbed 
phase. In 1 -|- 1 dimensions, this feature manifests in 
the fact that the length of the chains increases when de- 
creasing the variance of the probability distribution. In 
2-1-1 dimensions, we have reported the existence of a 
percolation transition which can be identified due to the 
presence of a peak in the susceptibility. In all cases, the 
jamming limit increases significantly when the probabil- 
ity distribution function becomes sharper, and is close, 
in the limit of a sharply peaked distribution, to the close- 
packing regime. 

Our model could caricature some real situations, not 
addressed by the sequential adsorption models proposed 
up to now, sharing in common the fact that the final po- 
sition of an adsorbed particle may exert an influence in 
the location of the next particle. Cases in which the par- 
ticles interact, when the kinetics exhibits memory effects 
or, in general, when the adsorption process may be in 
some way controlled by an external agent and therefore 
is not totally random, could accordingly find a descrip- 
tion in the model we have proposed. The results we have 




FIG. 8. Largest cluster of connected particles in 2 -|- 1 di- 
mensional correlated adsorption, a) u = oo; b) a = 2.8; c) 
a = 2.4, d) <T = 2.0. 



presented may provide new possibilities to the already 
existing ones about modelization of adsorption kinetics 
by sequential processes. 
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